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The standard procedure for detection of gravitational wave coalescing binaries sig- 
nals is based on Wiener filtering with an appropriate bank of template filters. This 
is the optimal procedure in the hypothesis of addictive Gaussian and stationary 
noise. We study the possibility of improving the detection efficiency with a class 
of adaptive spectral identification techniques, analyzing their effect in presence of 
non stationarities and undetected non linearities in the noise. 



1. Adaptive whitening filter and non stationarities 

We can define the output of an interferometric detector x(t) as the sum of 
noise n(t) and signal h(t, 0) of known shape and unknown parameters 0: 



It is well established that the optimal filter for the detection of the signal 
hit, 0), if the noise is Gaussian and stationary, is the Wiener matched filter 1 : 
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x(t) 



n{t)+h(t, 0). 
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being S(v) the noise Power Spectral density (PSD) and h(u, 9) the template 
for the signal we are looking for. 

If the noise is not stationary we have problems in estimating the PSD. 
We can avoid these by generalizing in time domain the Fourier space the 
"whitening" procedure in the Wiener filter, using an adaptive filter: with 
the "whitening" we transform a colored data sequence in a white one, i.e. 6- 
correlated with a "flat" PSD. We will use a linear filter produced by fitting 
the PSD with and Auto Regressive (AR) model 2 ' 3 ' 4 : 

p 

x[n] = a,kx[n — k] + <jw[n] (3) 
fe=i 

where ak,cr are the parameters of the model, P the order of the model and 
w[n] a driving normal white noise. The idea of adaptive filters is to find 
the optimal weights which minimize a cost function error e. For Least 
Squares based methods this is 

n 

e(n)=J2^e 2 m (4) 

i=l 

where e[i\n] is the error done in the prediction at time n with the at cal- 
culated at time n. The parameter A, where < A < 1, is called "forgetting 
factor" and is introduced to let the filter follow the changes in the pa- 
rameters of the noise model 4 . The adaptive filter we used, obtained by 
minimizing eq. 4, is the Least Squares Lattice (LSL) filter 4 . When we 
adaptively identified the parameters and their time variation with the LSL 
filter, we can use this filter to whiten the data: x[n] becomes the input 
colored sequence which enters in the filter, and w[n] is the output whitened 
sequence. We simulated a Virgo-like 5 non stationary noise by fitting the 
afe,CT parameters on the theoretical PSD for thermal and shot noise of the 
Virgo interferometer 2 and adding a sinusoidal modulation to a: 

cr[n] = a(l + Asm2irvn) (5) 

With the forgetting factor A less than one we succeeded in following the 
non-stationarity and in whitening the non stationary process obtaining a 
flat PSD as it is expected for a white process. In the first plot of Fig. 1 
we compare the standard deviation of the simulated non stationary time 
series with the standard deviation estimated by LSL filter, whereas in the 
second plot the PSD of the simulated non stationary time series and the 
PSD of the whitened time series are shown. The use of the adaptive filter 
with A < 1 is important in view of the occurrence of unknown slow non 
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stationarity in the interferometer data, where we cannot apply anymore the 
optimal filter. 



Figure 1. Top plot: standard deviation a of the non-stationary data -A = 0.5 and 
v = 0.1- and adaptive estimated standard deviation: the LSL filter is able to follow the 
non stationarity with A = 0.997. Bottom plot: PSD of the non stationary data before 
and after whitening. 



2. Performance of whitening and non-linearities 

Apparent non stationarities can appear as the effect of a stationary non 
linear process. For example, if we observe a spike, this can be a signature 
of both non stationarity or non linearity, whereas a linear stationary process 
has a low probability to be "spiky" . The whitening filter should then show 
a good performance also in presence of non linearities, because we may not 
be able to know which type of process we are analyzing. 

We simulated a non linear process for a Virgo-like model with an AR 
model with conditional heteroscedasticity (ARCH) 6 . The time evolution of 
a Virgo- like ARCH model is reproduced in Fig. 2, first plot: it presents a lot 
of spikes, but its PSD (on the right) is almost the same of the linear Virgo- 
like AR model, as the PSD retains only the linear features of a process. 
The whitening filter with A = 1, when applied to this process (plots on 
the second line of Fig. 2) does not completely eliminate the spikes, as the 
non-linearities are averaged away in the limit of an infinite number of data 
and are therefore irrelevant in the spectral estimation: the PSD is flat and 
the process has successfully been whitened, but the spikes could give rise 
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Figure 2. Performance of the adaptive filter in presence of non linearities. The plots on 
the right are the time series obtained with: the ARCH model, the same model filtered 
with A = 1 and with A = 0.987. The plots on the right are the PSD in the three different 
cases. 

to false alarms in the detection process, being identified as signals. For 
A = 0.987 (Fig. 2, third line) we find a reduction of the spikes. This can 
help in reducing the false alarm probability. 



3. Conclusions 

We showed the important role of forgetting factor in the adaptive LSL 
filter in identifying non stationarities of time series processes. Moreover 
this adaptive whitening procedure is apparently effective in reducing both 
non stationarity and non linearity effects. This could be important in order 
to improve the robustness of algorithm for gravitational wave detection 7 . 
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